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Resumen 


Buendía-Espinoza, J. C., Pérez-Rodríguez, P., González- 
Camacho, J. M., Pérez-Elizalde, S., Exebio-García, A., 
Rosengaus-Moshinsky, M., éz Peña, M. (julio-agosto, 2017). 
Identificación de cambios en la ciclogénesis del Atlántico 
Norte mediante un modelo de mezclas Gaussianas. Tecnología 
y Ciencias del Agua, 8(4), 5-18. 


Varios modelos climáticos sugieren que la frecuencia e 
intensidad de los ciclones tropicales cambiarán a finales 
del siglo XXL, como consecuencia del calentamiento global. 
Existen diferentes reportes científicos sobre el calentamiento 
global y su relación con la intensidad de los ciclones 
tropicales. Sin embargo, poco se ha investigado sobre el 
impacto del calentamiento global sobre la ciclogénesis en 
las diferentes cuencas oceánicas. En este trabajo se estima 
el número de regiones ciclogénicas y sus centroides en la 
cuenca oceánica del Atlántico Norte entre los periodos 1951- 
1975 versus 1976-2013 y 1951-1989 versus 1990-2013 mediante 
un modelo estadístico de mezclas finitas Gaussianas, para 
identificar si existen cambios significativos debidos al 
cambio climático. En el presente estudio, los cambios de 
un intervalo a otro se le atribuirán al cambio climático. La 
estimación de los parámetros de la función de densidad de 
probabilidades (fdp) de los componentes de las mezclas 
se hizo a través del algoritmo esperanza-maximización 
(EM). Las fdp se compararon mediante la distancia de 
Bhattacharyya y se estimó el percentil 95 a través de la 
técnica de remuestreo paramétrico. Los resultados obtenidos 
muestran que solamente hay dos regiones ciclogénicas en 
los dos intervalos de estudio, es decir, no hubo cambio en 
el número de regiones de un periodo a otro. Un segundo 
resultado es que existen diferencias estadísticamente 
significativas en la ubicación de los centroides de las regiones 
de génesis estimadas, sugiriendo un impacto debido al 
cambio climático de acuerdo con los datos analizados. 


Palabras clave: modelos de mezclas, verosimilitud, regiones 
ciclogénicas. 


Abstract 


Buendía-Espinoza, J. C., Pérez-Rodríquez, P., González- 
Camacho, J. M., Pérez-Elizalde, S., Exebio-García, A., Rosengaus- 
Moshinsky, M., € Peña, M. (July-August, 2017). Identification 
of changes in the North Atlantic cyclogenesis using a Gaussian 
mixture model. Water Technology and Sciences (in Spanish), 
8(4), 5-18. 


Several climate models suggest that the frequency and intensity 
of tropical cyclones will change at the end of the 21st century, as 
consequence of global warming. Several scientific reports have 
described how global warming affects the intensity of tropical 
cyclones. However, little research has been done on the impact of 
global warming on cyclogenesis in different ocean basins. In this 
work, the number of cyclogenetic regions and their corresponding 
centroids are estimated in the North Atlantic Ocean basin for the 
intervals 1951-1975 versus 1976-2013 and 1951-1989 versus 
1990-2013 through a finite Gaussian mixture model, to determine 
whether significant changes exist. In this study, the change from 
one interval to another will be attributed to climate change. The 
parameter estimation of the probability density function (fdp) of the 
mixture components was done using the Expectation-Maximization 
algorithm (EM). The fdp were compared using the Bhattacharyya'“s 
distance and the 95' percentile was estimated using the technique 
of parametric bootstrap. The results show that there are only two 
cyclogenic regions in each of the intervals that is no increase or 
decrease in the number of regions. A second result indicates that 
there are significant differences in the centroid's locations from one 
interval to another, suggestion an impact due to climate change 
according to the analyzed data. 
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Introducción 


Existe un interés creciente en conocer el im- 
pacto del cambio climático en la actividad de 
los ciclones tropicales, pues son unos de los 
fenómenos climáticos más destructivos entre 
los fenómenos naturales que provocan mayores 
desastres. A medida que la temperatura global 
promedio de la superficie del planeta se incre- 
menta, se espera que la intensidad, frecuencia, 
trayectorias, ubicación de ocurrencia y llegada a 
tierra de estos fenómenos se alteren por el clima 
actual. Walsh (2004) menciona que aunque no 
hay en este momento cambios perceptibles en 
las características de los ciclones tropicales que 
razonablemente podrían atribuirse al calenta- 
miento global, las predicciones de los modelos 
de circulación general muestran incrementos 
en su intensidad máxima entre 5 y 10% para 
el año 2050. También indica que las regiones 
de formación de huracanes probablemente no 
cambien, y que ha habido poco consenso respec- 
to a la modificación en el número de ciclones o 
trayectorias, además de la incertidumbre en las 
predicciones climáticas por algunas deficiencias 
en los modelos de circulación general; por lo 
tanto, si las predicciones de intensidades son co- 
rrectas, sus cambios se detectarán en el Atlántico 
después de 2050. Con base en los registros de las 
mejores trayectorias del Centro de Advertencia 
de Tifones de la Marina Estadounidense (Joint 
Typhoon Warning Center, JTWC) y el Centro 
Nacional de Huracanes (National Hurricane 
Center, NHC) de la Agencia Atmosférica y 
Oceanográfica Nacional (National Oceanogra- 
phic and Atmospheric Administration, NOAA), 
Emanuel (2005), Webster, Holland, Curry y 
Chang (2005), y más recientemente Elsner, 
Kossin y Jagger (2008) demostraron que la 
intensidad histórica de las tormentas ha aumen- 
tado tanto en el Pacífico Noroccidental (WNP) 
como en el Atlántico Norte. Webster et al. (2005) 
observaron un aumento considerable en los ci- 
clones tropicales en todas las cuencas oceánicas 
durante los últimos 30 años para las categorías 
más fuertes (4 y 5 de acuerdo con la escala de 
Saffir-Simpson); en particular, en las regiones 
ciclogénicas del Pacífico Noroccidental y del 


+» ISSN 0187-8336 


Atlántico Norte; 25 y 20%, respectivamente, de 
estas tormentas se presentaron en el intervalo 
de 1975-1989, y 41 y 25%, respectivamente, en 
el intervalo 1990-2004; esto significa que hubo 
aumentos de 16 y 5%, respectivamente. Sin em- 
bargo, los resultados obtenidos para la cuenca 
oceánica del Pacífico Noroccidental han sido 
cuestionados, pues la aparente tendencia que 
se Observa en la intensidad de los ciclones es 
parte de un gran oscilación interdecadal (Chan, 
2006) o de posibles errores de medición en el 
conjunto de datos (Knaff €: Zehr, 2007). Además, 
Klotzbach (2006), con base en el análisis de los 
registros de las mejores trayectorias del periodo 
1986-2006, encontró que la tendencia de la inten- 
sidad de los ciclones tropicales para la cuenca 
del Atlántico Norte no muestra evidencia de que 
haya cambiado, y la tendencia para el Pacífico 
Noroccidental presenta una baja considerable. 
También Kossin, Knapp, Vimont y Harper 
(2007), con base en el análisis de los registros de 
las mejores trayectorias, encontraron que no hay 
un aumento en la intensidad de los huracanes 
en cualquier cuenca distinta a la del Atlántico 
Norte en las dos décadas (1985-2005). 

Respecto a la ubicación de ocurrencia (géne- 
sis) y la lisis de los ciclones tropicales, también 
han experimentado variaciones. Knutson et al. 
(2010) mencionan que no hay ninguna prueba 
concluyente de que los cambios observados en 
la génesis, trayectorias, duración e inundaciones 
de oleaje de los ciclones tropicales superen la va- 
riabilidad esperada a partir de causas naturales 
en los últimos 50 años; sin embargo, mencionan 
que hay sugerencias de cambios en la ubicación 
de génesis y en las trayectorias de las tormentas 
observadas en el océano Atlántico ofrecidas 
como una explicación de la falta de tendencias 
crecientes de la llegada a tierra de ciclones en los 
Estados Unidos y la Costa del Golfo. Asimismo, 
Kossin, Emanuel y Vecchi (2014) encontraron 
que los cambios observados en el cizallamiento 
vertical del viento y la intensidad potencial 
proporcionan evidencias o pruebas de que el 
desplazamiento global de los ciclones tropicales 
hacia los polos en ambos hemisferios, es decir, 
fuera de los trópicos, está siendo modificado por 
cambios ambientales sistemáticos en los últimos 
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30 años. También mencionan que si estos cam- 
bios ambientales continúan, un desplazamiento 
simultáneo concomitante hacia los polos de la 
latitud donde los ciclones tropicales alcanzan 
su intensidad máxima tendría consecuencias 
potencialmente profundas sobre los residentes 
costeros e infraestructura. Los ciclones tropica- 
les también tienen un papel importante en el 
sostenimiento de los recursos hídricos regiona- 
les; de igual forma, un desplazamiento hacia 
los polos de las tormentas podría amenazar 
el abastecimiento de agua potable en algunas 
regiones, al tiempo que aumentan los eventos 
de inundación, entre otras consecuencias. Por su 
parte, Mori, Kuniyoshi, Nakajo, Yasuda y Mase 
(2013) simularon con un modelo estadístico el 
impacto del calentamiento global sobre los cen- 
troides de la ciclogénesis y de la ciclolisis de las 
diferentes regiones ciclogenéticas. Concluyeron 
que éstos se desplazarán hacia el centro de las 
cuencas oceánicas y que los cambios futuros en 
las condiciones dinámicas y termodinámicas 
en los océanos influirán en la frecuencia de la 
génesis de los ciclones tropicales y en el cambio 
de localización hacia el centro de las diferentes 
cuencas. Cabe mencionar que las temperaturas 
de la superficie del mar son más cálidas en 
el centro del océano que en las orillas (Chan, 
2007; Yokoi $: Takayabu, 2009), y que un ciclón 
que se desarrolla en la parte central dura más 
tiempo. Bajo tal contexto, estos autores men- 
cionan que los cambios en la intensidad de los 
ciclones tropicales están más relacionados con 
el desplazamiento de los centroides que con la 
modificación de la temperatura de la superficie 
del mar. Holland y Bruyere (2013) desarrollaron 
un Índice de Cambio Climático Antropogénico 
(ACCI, por sus siglas en inglés, Anthropogenic 
Climate Change Index) para investigar la con- 
tribución del calentamiento global a la activi- 
dad de los ciclones tropicales. Dichos autores 
encontraron que en todas la cuencas oceánicas 
la frecuencia anual de los ciclones o huracanes 
tropicales no se ve afectada por el cambio antro- 
pogénico. Sin embargo, los huracanes categoría 
4 y 5 han aumentado a tasas de 25 a 30% por 
grado Celsius ("C) de calentamiento global; en 


contraste, los huracanes de categorías 1 y 2 han 
disminuido a tasas similares. 

Además de la variación en la intensidad y 
la génesis, las precipitaciones inducidas por los 
ciclones tropicales también han experimentado 
cambios a largo plazo con la oscilación inter- 
decadal. Kim, Ho, Lee, Jeong y Chen (2006), y 
Lau, Zhou y Wu (2008) encontraron un cambio 
sustancial en la precipitación durante la estación 
de ciclones, manifestado como un aumento en 
eventos y cantidad de lluvia extrema tanto en 
el Pacífico Noroccidental como en el Atlántico 
Norte. Sin embargo, en lo que respecta a las re- 
laciones entre los ciclones tropicales y las lluvias 
extremas, descubrieron diferencias centrales 
entre los dos dominios oceánicos: en el Atlántico 
Norte, las contribuciones de los ciclones tropica- 
les a la precipitación extrema son mayores que 
en el Pacífico Noroccidental. 

En este trabajo de investigación se aplicó 
un modelo estadístico de mezclas Gaussianas 
(sección 2). Por una parte, para determinar el 
número de regiones ciclogénicas; por otra, para 
determinar los cambios temporales y espaciales 
de los centroides de dichas regiones en la cuenca 
oceánica del Atlántico Norte para dos periodos 
(1951-1975 versus 1976-2013 y 1951-1989 versus 
1990-2013). Se utilizaron los puntos de inicio 
de los datos (latitud y longitud) de las “mejores 
trayectorias” o IBTrACS (Ncdc.noaa.gov, 2017). 
Se determinaron las funciones de densidad de 
probabilidades (fdp) de las regiones de génesis 
a través del algoritmo esperanza-maximización 
(EM), las cuales se evaluaron y compararon para 
verificar los cambios espacio-temporales. En la 
sección 3 se muestran los resultados, y su rela- 
ción e interpretación con el fenómeno natural en 
estudio. Finalmente, en la sección 4, se concluye 
con una breve discusión. 


Materiales y métodos 


Descripción y alcance geográfico de la base de 
datos de ciclones tropicales 


El análisis de los datos fue de 1951-2013 (ver 
figura 1), dividido en dos periodos: 1951-1975 
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Figura 1. Distribución espacial de los puntos de ubicación de ocurrencia de los ciclones tropicales en la región del Atlántico 


Norte desde el año 1951 hasta 2013. Fuente: elaboración propia con base en los datos del IBTrACS. 


versus 1976-2013 y 1951-1989 versus 1990-2013. 
El primero se determinó con base en el calen- 
tamiento ocurrido durante el siglo XX en los 
intervalos de 1928-1975 y 1976-2013 (ver figura 
2). Cabe mencionar que se tomó como valor 
de inicio 1950, porque a partir de este año fue 
usado de manera rutinaria el reconocimiento 
aéreo para supervisar los ciclones tropicales. 
Esto significa que la información sobre ciclones 
medidos antes de esta fecha es menos confiable 
(Vecchi £: Knutson, 2008, citados por Villarini, 
Vecchi, éz Smith, 2011). 

Respecto al segundo periodo, se determinó 
de acuerdo con Webster et al. (2005), quienes 
observaron un aumento considerable en los 
ciclones tropicales en todas las cuencas oceá- 
nicas durante los últimos 30 años para las 
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categorías más fuertes (4 y 5 con base en la 
escala de Saffir-Simpson) debido al incremento 
de la temperatura de la superficie marina (ISM). 
Cabe mencionar que dicho intervalo de tiempo 
lo dividieron en dos periodos: 1975-1989 versus 
1990-2004. 


Modelo de mezclas Gaussianas (MMG) 


Bajo el modelo de mezclas finitas de distribu- 
ciones normales, ajustado en este estudio, cada 
punto de ubicación de ocurrencia x, = (latitud, 
longitud ) de los ciclones tropicales se consi- 
dera que proviene de una súper población, la 
cual es una mezcla de un número finito K de 
poblaciones con proporciones O pesos K,,... T,, 
respectivamente, donde: 
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a) Variaciones de la temperatura media mundial de la 
superficie del mar para el intervalo 1870-2008. 
Fuente: Knudsen et al. (2001) 
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Fuente: Knudsen et al. (2001) 


Figura 2. Variaciones de la temperatura media de la superficie del mar global y del Atlántico Norte. 


> a=l; n.=0 (1) 


La función de densidad probabilidad (fdp) 
de la variable aleatoria X (de dimensión d) en la 
forma de mezclas finitas es: 


p(x,0)= xy p (150) 
E AS) 
2% 2 o 2) 


donde q denota el vector de parámetros de la 
mezcla; p(x; 0) es la función de densidad de 
probabilidades correspondientes al K-ésimo 
componente; x, es el peso de la mezcla; y 0, 
denota el vector de todos los parámetros desco- 
nocidos de los elementos de los vectores de las 
medias u, y elementos distintos de matrices de 
covarianzas 2, para¡=1,..., K, que pertenecen 
a algún espacio de parámetros 6. 


a) Estimación de los parámetros utilizando el 
algoritmo EM 


Existen varios procedimientos para determi- 
nar los parámetros de un modelo de mezclas 
Gaussianas (MMG) de un conjunto de datos 
(McLachlan éz Basford, 1988). Sin embargo, el 


método más popular es el de máxima verosi- 
militud (ML, por sus siglas en inglés). En este 
trabajo, la estimación de los parámetros se hizo 
utilizando técnicas de máxima verosimilitud 
vía el algoritmo EM (Dempster, Laird, € Ru- 
bin, 1977; Redner € Walker, 1984; Maitra ézx 
Melnykov, 2010). 

El algoritmo EM para mezclas Gaussianas es 
un proceso iterativo que consiste de dos etapas: 
a) esperanza (E) y b) maximización (M). En el 
caso de componentes Gaussianos, la densidad 
de la mezcla contiene los siguientes parámetros: 
TL, U, y >, donde¡=1,..., K. La log-verosimilitud 
condicionada completa esperada para el conjun- 
to de datos, conocida como función Q, es igual a: 


Q(9,%,,--,x,) = 


1HÉ Tod 
37 (In]2/+(x,-1) El (1) 


¡A j=1 


¿SY im(x)- 2 im(25) (8) 


2 j=1 


El paso-E consiste en la actualización de las 
probabilidades de pertenencia x,'*”, dada la 
estimación de los parámetros actuales d”: 


li > ap (x aa (4) 
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El paso-M tiene una solución analítica: 


n + + + d 
se, Sa 6 ie Mit ul ») d 
j Na 
ia Y 


El algoritmo EM comienza con una esti- 


mación inicial del valor de los parámetros q, 
llamada q”. Luego, mediante las ecuaciones 
(4) y (5), se estiman los valores de los nuevos 
parámetros, llamado q'*”. El proceso se repite 
hasta que la diferencia entre dos evaluaciones 
sucesivas de la log-verosimilitud sea menor que 
una épsilon dada e, es decir: 


IS 


en este trabajo se fijó e = 1 x 107. Este resultado 


<E (6) 


depende de la selección de los parámetros ini- 
ciales (Seidel, Mosler, éz Alker, 2000). 


b) Inicialización en el algoritmo EM 


El algoritmo EM es un procedimiento iterativo 
de maximización que depende del valor inicial 
de los parámetros, pues la función de verosimi- 
litud puede tener máximos locales (McLachlan 
éz Peel, 2000). Por lo tanto, una buena iniciali- 
zación es crucial para encontrar los estimadores 
de máxima verosimilitud. 

Se han sugerido diferentes procedimientos 
de inicialización en la literatura (Figueiredo éz 
Jain, 2000; Maitra, 2009); sin embargo, ningún 
método supera a los demás. En este trabajo 
se utilizó el procedimiento de Fraley, Raftery, 
Murphy y Scrucca (2012), implementado en la 
biblioteca de funciones de R (R Core Team, 2016) 
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mclust, para encontrar los valores iniciales de 
los parámetros que permiten obtener el valor 
máximo en el marco de mezclas Gaussianas 
multivariadas. 


Identificación del número óptimo de 
componentes o grupos 


Hay una vasta lista de literatura dedicada al 
tema de la elección de K (número de compo- 
nentes). McLachlan y Peel (2000) proporcionan 
una interpretación detallada de los diferentes 
enfoques disponibles para abordar este pro- 
blema. La mayoría de los métodos destinados 
a la estimación de K se divide por lo general 
en dos categorías: modelos basados en el prin- 
cipio de la parsimonia y modelos basados en 
procedimientos de prueba, ambos sustentados 
en la función de log-verosimilitud. Sin embargo, 
en este estudio, K se determinó mediante un 
método heurístico, conocido como partición 
alrededor de los medoides (PAM, por sus siglas 
en inglés, Partitioning Around Medoids). 

El algoritmo de la PAM se basa en la forma- 
ción de K particiones u objetos representativos 
(medoides) de n observaciones de un conjunto 
de datos. Un medoide se define como la obser- 
vación de un agrupamiento, cuya diferencia pro- 
medio, con respecto a todas las observaciones en 
el grupo, es mínima. Se eligen aleatoriamente K 
medoides de un conjunto de datos. El medoide, 
que representa un grupo, se ubica en el centro 
del grupo. Los objetos restantes se agrupan con 
el medoide al que son más similares, basándo- 
se en la distancia entre el objeto y el medoide. 
La estrategia, entonces, es reemplazar uno de 
los medoides por los no medoides, siempre y 
cuando la calidad del agrupamiento mejore. 
Esta calidad es estimada usando una función 
de costo que mide el promedio de disimilaridad 
o diferencia entre un objeto, y el medoide de su 
grupo (Kaufman é: Rousseeuw, 2005). 

El método PAM genera una gráfica, conocida 
como gráfica de “siluetas”. Para cada obser- 
vación se muestra una medida que indica la 
calidad de la clasificación. Valores cercanos a 1 
indican que la observación está bien clasificada 
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en su grupo; valores cercanos a 0 significan que 
la observación podría pertenecer a otro grupo; 
en valores cercanos a -1, la observación está 
pobremente clasificada. En la gráfica hay una 
medida resumen denominada ancho promedio 
de silueta o coeficiente de silueta, que se inter- 
preta de acuerdo con el cuadro 1. 


Medida de la distancia entre MMG 


De acuerdo con Sfikas, Constantinopoulos, 
Likas y Galatsanos (2005), y Fukunaga (1990), 
la distancia de Bhattacharyya puede ser utiliza- 
da para comparar fdp de modelos de mezclas 
Gaussianas y además tiene una expresión en 
forma cerrada. Esta distancia se utilizó para me- 
dir la distancia entre los grupos de los ciclones 
tropicales para los dos intervalos de estudio y 
se define como sigue: 


alo E J (u,-1,) 


1 2 
—In E 
a (7) 


donde Uy 2,y U, Za corresponden a los vecto- 
res de medias, y las matrices de varianzas y 
covarianzas de los núcleos de las densidades 
Gaussianas, respectivamente. 


Estadístico de prueba para la comparación de 
las fdp del MMG 


Se aplicó la técnica paramétrica de remuestreo 
para obtener el percentil 95 de la distribución 


empírica del estadístico de prueba, que es la dis- 
tancia de Bhattacharyya entre los componentes 
de los modelos de mezclas Gaussianas, para 
verificar si existen diferencias estadísticamente 
significativas entre ellos (Engel, 2010). 

El método de remuestreo (Efron, 1979; Efron 
£z Tibshirani, 1993) se usa para encontrar inter- 
valos de confianza en situaciones donde es im- 
posible obtener de manera analítica la distribu- 
ción muestral del estimador (para mayor detalle, 
ver DiCiccio €: Tibshirani, 1987; Hall, 1988). Esta 
técnica se sustenta teóricamente en dos conside- 
raciones: 1) la función de distribución verdadera 
F(x) se estima mediante la función de distribu- 
ción empírica É(x) de acuerdo con el tegrema 
Glivenko-Cantelli, que muestra que F(x) > F(x) 
conforme 1 => co (Bickel éz Freedman, 1981); y 2) 
de acuerdo con la propiedad de consistencia, la 
F(0) de una muestra dada puede aproximarse 
mediante la distribución muestral del remues- 
treo É'(6') cuando el número de remuestreos 
es grande y puede también aproximarse Elx) a 
F(x). Bajo estos supuestos, Babu y Singh (1983) 
demuestran que FS (8) ÉE(0), cuando el número 
de remuestreos es suficientemente grande. 


Procedimiento de remuestreo 


1. Se genera una variable aleatoria U-Uniforme 
(0,1). 

2. Si ve[Y/,1,),. a) donde 1, 
corresponde a la probabilidad del ¡-ésimo 
componente del modelo de mezclas, 
entonces se generan variables aleatorias 
a partir de la distribución del ¡-ésimo 
componente. Dicho componente tiene una 
distribución normal bivariada, obtenida 
mediante el algoritmo EM. 


Cuadro 1. Rango del coeficiente en silueta. 


Rango del coeficiente de silueta 


Interpretación 


0.71-1.00 Estructura fuerte 

0.51-0.70 Estructura razonable 

0.26-0.50 Estructura débil y podría ser artificial 
<0.25 Ninguna estructura sustancial 


Fuente: Kaufman y Rousseeuw (2005). 


ISSN 0187-8336 + 


Tecnología y Ciencias del Agua, vol. VII, núm. 4, julio-agosto de 2017, pp. 5-18 


Tecnología y Ciencias del Agua, vol. VII, núm. 4, julio-agosto de 2017, pp. 5-18 


Buendía-Espinoza et al., Identificación de cambios en la ciclogénesis del Atlántico Norte mediante un modelo de mezclas Gaussianas 


3. Serepiten los pasos 1 y 2 hasta que se tenga 
la cantidad deseada de muestras de la 
mezcla de la distribución. 

4. Se definen las variables aleatorias X, — N, (u, 
Y), ¡=1, ..., n obtenidas en el paso anterior 
como muestra aleatoria. 

5. Se obtiene una muestra bootstrap X' = (x,, 
..., X,) muestreando aleatoriamente con 
reemplazo n veces los datos originales x,, 
..., X,; el tamaño de la muestra aleatoria es el 
mismo que el de la muestra de remuestreo, 
y las X tienen probabilidad n”, siendo igual 
en cada una de las X.. 

6. Se calcula el estadístico d, (f,2) (distancia 
de Bhattacharyya) de este remuestreo, 
produciendo d, (f,3). 

7. Se repite el paso 4 B veces. La ley de los 
números grandes indica que si B es lo 
suficientemente grande, se obtendrá una 
buena aproximación a la densidad verda- 
dera d,(f, 3). En este trabajo se fijó B=1 000. 

8. Se construye la distribución de probabilidad 
de las Bd; (f,g), asignando probabilidad B” 
a cada d,, (f,2). Esta es la estimación de la 


distribución d, (f,g), E” (4 10 g)) 
Prueba de hipótesis a partir del remuestreo 


La técnica de remuestreo permite realizar la 
prueba de hipótesis de similitud de las fdp de 
los grupos de los modelos de mezclas Gaussia- 
nas entre los diferentes intervalos, es decir, H.,: 
d(f, 3) =0. La regla de decisión es rechazar H, si 
d,(f, 3) es grande. El procedimiento consiste en: 


1. De los datos originales de la muestra se 
obtienen los estimadores (1, 8) de los 
componentes del modelo de mezclas normal 
bivariado mediante el algoritmo EM. 

2. Mediante el procedimiento de remuestreo 
descrito anteriormente se obtiene d,(f,2) 
bajo H, (es decir, d, (f,9)y,,), realizando 1 000 
muestras de remuestreo de la distribución 
normal bivariada bajo H,; d,(f, g) =0. 

3. Cada remuestreo deberá ser del mismo 
tamaño que el de la muestra inicial. 
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4. Se calcula el estadístico d, (f,2) para 
cada muestra de remuestreo, y con ellos 
se construye la función de distribución 
empírica de d, (f,2). 

5. La prueba de hipótesis de la distancia de las 
fdp de los grupos es H.; d, (f, 2) con a. =0.05; 
esto equivale a obtener el percentil 95, y 
rechazar H, si d a (f,8) es mayor que éste. 


Resultados 


Estimación de la función de densidad de 
probabilidad 


De acuerdo con la figura 3, el número de gru- 
pos por intervalo (1951-1975 versus 1976-2013 y 
1951-1989 versus 1990-2013) es dos (es decir, K = 
2 componentes). Su ancho promedio de silueta 
fue 0.585 para ambos intervalos (ver figura 
4), lo cual significa que tienen una estructura 
razonable (ver cuadro 1). 

Los valores de los parámetros iniciales para 
cada uno de los componentes de las mezclas 
(K = 2) se determinaron mediante el algoritmo 
implementado por Fraley et al. (2012). El vector 
de parámetros q se estimó iterativamente me- 
diante el algoritmo EM. La función de densidad 
estimada para cada uno de los modelos de 
mezclas Gaussianas se ajusta a la distribución 
espacial de los puntos de ubicación de ocurren- 
cia de los ciclones tropicales en ambos intervalos 
(figura 5). Los puntos que están muy alejados 
de los centroides exhiben un comportamiento 
atípico desde el punto de vista estadístico, pero 
dado que se sabe que son datos reales, se de- 
cidió mantenerlos como muestras legítimas al 
momento de realizar los análisis. Otra opción 
para manejo de este tipo de datos es el uso de 
mezclas con colas pesadas, como por ejemplo la 
distribución t-bivariada. 

Además, en la figura 6 se muestra la ubica- 
ción de los centroides de cada uno de los grupos 
para ambos periodos de estudio, así como su 
función de densidad estimada por periodo. Los 
puntos negros corresponden a los datos. Los dos 
grupos en los dos periodos de estudio aparente- 
mente indican que tienen la misma orientación. 


Buendía-Espinoza et al., Identificación de cambios en la ciclogénesis del Atlántico Norte mediante un modelo de mezclas Gaussianas 


mn 
[« LO 
< 15 
== 0 
g 
¡0 
3 
E 2 
3 S 
o 
io] 
8 
T Y 
E [=] 
2] 
E 
a 
9 =] 
DN 
[= o 
< 
5 10 15 20 
Núm. de grupos 
(a) Intervalo 1951-1975 
E 
< 8 
a (= 
g 
¡0 
3 
M5] 
< 3 
2 a 
2 
us] 
E a 
K 
g Ss 
a, 
e] 
al 
(a) 
SS 
¿ $ 
S 


Núm. de grupos 


(c) Intervalo 1976-2013 


0.45 0.50 0.55 0.60 


Ancho promedio de la silueta (APS) 


0.40 


5 10 15 20 


Núm. de grupos 


(b) Intervalo 1951-1989 


0.45 0.50 0.55 


Ancho promedio de la silueta (APS) 


0.40 


Núm. de grupos 


(d) Intervalo 1990-2013 


Figura 3. La gráfica número de grupos versus amplitud de silueta. Esta gráfica muestra que a medida que aumenta el número 


de grupos por intervalo disminuye el APS en las cuatro gráficas, pero disminuye la consistencia en su estructura de acuerdo 


con el cuadro 1. 


Prueba de remuestreo no paramétrica para 
comparar las fdp de los MMG 


La obtención del percentil 95 de la distribu- 
ción empírica del estadístico d, (f,3) equivale 
a rechazar H, si d, (£,2) > 9.95: Por lo que con a 
= 0.05 y d, (f,g) = 0.081 se rechazó la hipótesis 
nula, indicando que existen diferencias esta- 
dísticamente significativas entre las fdp de las 
mezclas. 


Discusión y conclusiones 


Con base en los resultados obtenidos, en cada 
uno de los periodos de estudio se encontraron 
dos regiones de génesis de ciclones tropicales 
en la cuenca oceánica del Atlántico Norte; esto 
significa que no existen indicios de que se es- 
tén formando nuevas regiones ciclogenénicas 
debido al cambio climático. Por lo tanto, se 
concluye que solamente hay dos regiones que 
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n=341 2 grupos C; 


1: 19510.64 


2: 146 10.50 


0.0 0.2 0.4 0.6 0.8 1.0 
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Ancho promedio de la silueta: 0.58 


(a) Intervalo 1951-1975 


n = 585 2 grupos Cj 


1: 33310.60 


2: 252 10.58 


0.0 0.2 0.4 0.6 0.8 1.0 


Ancho de la silueta, S; 
Ancho promedio de la silueta: 0.59 


(c) Intervalo 1976-2013 


j: n¡| promedio iec;Sj 


j: nl promedio iec;9 ; 


n=560 2 grupos Cj 


jon; | promedio iec;9j 


1: 32610.64 
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Ancho promedio de la silueta: 0.60 


(b) Intervalo 1951-1989 
n=366 2 grupos Cj 
j: 1; | promedio tec 


1: 20110.57 


2: 165 10.57 


AN: al ET, E a | 
0.0 0.2 0.4 0.6 0.8 1.0 


Ancho de la silueta, S; 


Ancho promedio de la silueta: 0.57 


(b) Intervalo 1990-2013 


Figura 4. El Ancho promedio de la silueta, de manera general, indica que hay una buena estructura en los grupos elegidos, con 


la mayoría de las observaciones que parecen pertenecer a la agrupación en que están de acuerdo con el cuadro 1. 


generan los ciclones tropicales desde 1951 
hasta el 2013. 

Los centroides de las dos regiones cicloge- 
nénicas ubicadas en la cuenca oceánica del At- 
lántico Norte han experimentado cambios en su 
localización en el intervalo de estudio; es decir, 
se han movido hacia el centro y hacia los polos 
de la cuenca, lo cual coincide con Knutson et al. 
(2010) y Kossin et al. (2014), quienes mencionan 
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que los cambios ambientales sistemáticos en los 
últimos 50 años han provocado que los ciclones 
tropicales se desplacen hacia los polos en am- 
bos hemisferios, esto es, fuera de los trópicos. 
Dichos desplazamientos encontrados en este 
trabajo también coinciden con las proyeccio- 
nes estadísticas hechas por Mori et al. (2013) 
mediante datos generados a través de modelos 
de circulación general para finales del siglo XXI 
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Figura 5. Función de densidad estimada por intervalo; a) muestra la marginal de la longitud y la latitud para los intervalos 
1951-1975 y 1951-1989; b) indica la marginal de la longitud y la latitud, pero para los intervalos 1976-2013 y 1990-2013. Los 


puntos negros son los puntos de los datos. Los dos grupos en los dos periodos de estudio aparentemente indican que tienen la 


misma orientación. 


Latitud 


= ».Periodo 1990-2013 


Latitud 
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Longitud Longitud 
(a) Periodo 1951-1975 vs. 1951-1989 (b) Periodo 1976-2013 vs. 1990-2013 


Figura 6. Regiones de génesis de los ciclones tropicales y sus centroides aparentemente se han movido hacia el centro 
y hacia los polos de la cuenca oceánica del Atlántico Norte en los dos intervalos de estudio. 
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respecto al impacto del cambio climático sobre 
el desplazamiento de las regiones de génesis 
de los ciclones tropicales. 

En específico, el desplazamiento del cen- 
troide de la región de génesis del Golfo de 
México del periodo 1951-1975 versus 1951-1989 
fue de 1.1942 grados al este y 0.1279 al norte 
(ver panel a) de la figura 6); mientras que 
el centroide de la misma región del periodo 
1976-2013 versus 1990-2013 se desplazó 3.0021 
grados al este y 0.3067 al sur (ver panel b) de 
la figura 6). Respecto al desplazamiento del 
centroide de la región del Atlántico-tropical 
del periodo 1951-1975 versus 1951-1989, éste 
fue de 1.25085 grados al este y 1.18242 al sur, 
en tanto que el desplazamiento del centroide 
de la misma región para el periodo 1975-2013 
versus 1990-2013 fue de 0.06608 grados al este 
y 1.15809 al sur. 

Además, el hecho de que los centroides de 
las regiones ciclogénicas se desplacen hacia el 
noreste y sur de la cuenca oceánica provocaría 
que la precipitación inducida por los ciclones 
tropicales sufriría cambios tanto en distribu- 
ción espacial como en eventos y cantidad, tal 
y como lo mencionan Kim et al. (2006) y Lau 
et al. (2008). Esto significa que en la región 
central del continente americano, zona donde 
se ubica México, su precipitación disminuiría 
ligeramente, como indican Houghton et al. 
(2001) en sus proyecciones de precipitación 
para el siglo XXI debidas al cambio climático. 

Chan (2007), y Yokoi y Takayabu (2009) 
indican que una consecuencia importante 
que se generaría por el desplazamiento de 
los centroides de la cuenca es que los ciclones 
tropicales durarán más tiempo debido a que 
se generarán más hacia el centro de la cuenca 
oceánica; las temperaturas de la superficie 
del mar en el centro son por lo general más 
cálidas que en las regiones externas. El alcance 
de nuestro estudio es limitado, ya que no es 
posible saber si el cambio en los puntos de 
ciclogénesis impactará sobre la duración de 
los mismos o los cambios en los puntos de 
llegada a tierra. 
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